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We present full quantum mechanical calculations of the hybridized plasmon modes of two 

'"^ nanowires at small separation, providing real space visualization of the modes in the transition 

i-H from the classical to the quantum tunneling regime. The plasmon modes are obtained as cer- 

4^ tain eigenfunctions of the dynamical dielectric function which is computed using time dependent 

(\j density functional theory (TDDFT). For freestanding wires, the energy of both surface and bulk 
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plasmon modes deviate from the classical result for low wire radii and high momentum transfer 
due to effects of electron spill-out, non-local response, and coupling to single-particle transitions. 



F^ For the wire dimer at separations 1.5-0.3 nm, the spatial shape and energy of the plasmon modes 

Q are continuously altered with distance and deviate from the classical result due to the increasing 
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overlap of unoccupied electronic states of the two wires (virtual tunneling). At separations around 
2-3 A electrons at the Fermi energy can tunnel between the two wires and a charge-transfer mode 
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I. INTRODUCTION chemical sensing, to photocatalysi^ and 

photovoltaics^. Applications of plasmons in 

Plasmons are collective oscillatory modes ^^^^ generation nano-electronics^ and quan- 

of the electron system that are found tum information technologjja have also been 

m solids, at extended surfaces, and at envisioned. Advances in nano-fabrication 

the surface of coniined metal structured, techniques and transmission electron mi- 

The localized surface plasmon resonances ^roscopy have recently enabled the study 

sustained by metalhc nanoparticles^ have ^f plasmon excitations on a subnanometer 

found widespread use in apphcations rang- ^^^g^^ ^^^^^ 9 j^ particular, the use of 
ing from molecular spectroscopy^ and 
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highly confined electron beams can provide of the plasmon modes of the individual sub- 
spatial resolution of the electron energy loss system^I'^Sl. Recent experiments on coupled 
spectrum (EELS) yielding unprecedented nanoparticles found the energy of hybridized 
information about the dynamical properties modes to be blue-shifted relative to the clas- 
of confined electron systems^^''^. sical result at separations less than 1 nmP^, 
When the size of a plasmonic structure and a semi-classical approach, using a quan- 
reaches the nano-scale, quantum effects be- tum corrected model, found that the onset of 
gin to play a dominant role and the widely quantum tunneling should occur at around 
used classical models of electronic response 0.3 nm particle separation^^. However, until 
become insufficient^^. In this size regime, now the inffuence of quantum effects on the 
the theoretical description must be refined to spatial shape of the hybridized modes has not 
account for quantum effects such as surface been explored. 

scattering, electron spill-out at the surface. In this work we apply a recently developed 
coupling of plasmons to single-particle exci- method^^ to compute and visualize the hy- 
tations, and the non-locality of the electronic bridization of the quantum plasmons of two 
response. We note that some of these effects, parallel metallic nanowires at separations be- 
in particular the non-local response and spill- tween and 2 nm. The quantum plasmon 
out, can be incorporated into semi-classical modes are defined as certain eigenfunctions of 
hydrodynamic model^^^^^*^, however, the va- the dynamical dielectric matrix, which we ob- 
lidity and predictive power of such an ap- tain using linear response TDDFT, and they 
proach in the true quantum regime is not represent the possible eigenmodes of the elec- 
fuUy clarified. tron system damped by coupling to single- 
One important hallmark of the quantum particle excitations, 
regime is electron tunneling which occurs We first show that the plasmon modes of 
when two systems are brought into close single wires can be classified as surface or 
proximity. In fact tunneling leads to a break- bulk modes with a specific angular symme- 
down of the local field enhancement predicted try index m and radial node index n. The 
by the classical theory at sub-nanometer maximum number of modes found for a given 
separations, and leads to the formation of wire radius is limited by Landau damping 
a charge-transfer plasmon .'^^'^. Classically, which increases with n and m. For the dou- 
the plasmon modes of a composite system ble wire system, we find hybridized plasmon 
can be obtained as even/odd combinations modes. For separations below L5 nm, the 



hybrid modes begin to change in shape and culated densities within the jeUium model, 
their energy deviates from the classical result. Electron spill-out extends approximately 3 A 
These effects are caused by changes in the outside the jellium edge, which is important 
energy and wave function of electronic states for a correct description of surface plasmon 
lying several electron volts above the Fermi modes in particular. The local density ap- 
level (virtual tunneling). At separations be- proximation (LDA) was used for exchange- 
low 2-3 A we enter the true tunneling regime correlation. The unit cell was sampled on 
where the single-particle states around the grid-points of 0.2 A spacing, and 20 k-points 
Fermi energy start to overlap leading to an were used to map the ID Brillouin zone along 
abrupt change in the energy and shape of the wire, except for the study small mo- 
the plasmon modes with certain modes dis- mentum transfers in the response calculation 
appearing and being replaced by a charge- where 80 k-points were used, 
transfer mode. The microscopic dielectric function, 

eg(r, r',CL;), of the system was calculated 

using; linear response TDDFT^ within the 

II. METHOD ^ ^ 

random phase approximation (RPA), using 

The DFT calculations have been car- a 15 eV energy cutoff for the reciprocal 

ried out with the electronic structure code lattice vectors perpendicular to the wire 

GPAW^SE3 using the Atomic Simulation and including electronic states 10 eV above 

EnvironmentP. We study metalhc Na wires the Fermi level. The dielectric function was 

with radii in the range 2 to 8 A modeled sampled on an energy grid ranging from 

within the jellium approximation. Here the to 10 eV with a spacing of 0.01 eV, and a 

ions are described with a constant positive broadening of 0.2 eV was apphed to smear 

background density that terminates abrupt the discreteness in energy. The wave-vector 

at the edges, while the electron density is along the wire, gy, was varied between 0.03 

converged self-consistently. Since the jelhum A and 0.8 A in order to study the energy 

approach corresponds to a true free-electron dispersion. Convergence have been checked 

metal, it is a good approximation to sim- with respect to aU parameters, 

pie metals such as Na. The wire was con- The plasmon eigenmodes were obtained 

structed within a periodic supercell with 28 by diagonalizing the dielectric function for 

A of vacuum perpendicular to the wire, see each point on the frequency gricP^. This 

Fig. [l} showing the unit cell and the cal- leads to a set of distinct eigenvalue- curves 
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FIG. 1. (a) Sketch of the unit ceh with gray region indicating the positive charge density of a 
cyhndrical jelhum wire, surrounded by vacuum 2v = 28 A. i? is varied from 1 to 8 A and the unit 
cell is 3 A along the wire (b) Cross-sectional view of the wire showing the fixed positive background 
density profile p+ (red) and the converged electron density p (blue), plotted relative to the free 
electron density of bulk Na, Phom- 



that evolves smoothly with frequency, each 
curve corresponding to a separate plasmon 
eigenmode, see Fig. [2] The corresponding 
eigenf unctions, ^mdl^^^n), defines the spa- 
tial form of the induced potential (and den- 
sity) of the eigenmodes. When the damp- 
ing of the plasmon modes is small (which 
leads to a vanishing imaginary part of the 
eigenvalues), the eigenmodes can be identi- 
fied at points where the real part crosses zero 
from negative to positive values. For larger 
damping, the maxima of the eigenvalue loss- 
spectrum, — Ime~-^(cL;, r, /), provide a more 
accurate measure of the energy of plasmon 
resonances. Since large damping is found 
in some cases, we have used the last ap- 
proach throughout, and taken the half width 
of half maximum of the peaks as a measure 
of the plasmon energy. The real part of the 



eigenvalue- curves each resembles a Drude- 
like behavior for the dielectric function, and 
has been fitted to a single-pole model: 



en{uj) = 1 



C^r 



(1) 



which is used to find the mode strength, a^- 
The parameter /5^ corresponds to the average 
energy of the single particle transitions con- 
tributing the the given eigenmode, and 7^ is 
the damping of the mode. 

III. RESULTS 

The obtained plasmon eigenmodes for sin- 
gle wires of radius 2, 5 and 7.5 A is shown 
in Fig. [3] (a)-(c), where the induced poten- 
tial (first row) and density (second row) is 
plotted for increasing energy going left to 
right. The modes are characterized as sur- 
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FIG. 2. (a) Real part of the eigenvalues of the 
dielectric function eq{(jj^r^r') for a 5 A radius 
jellium wire. Each eigenvalue-curve that extends 
to negative values defines the existence of a plas- 
mon eigenmode, and the crossing with zero from 
below (red points), approximates the plasmon 
energy of the current mode, (b) For a more 
correct description of the mode energies in the 
presence of large damping, maxima in the eigen- 
values of — Ime~-'^(r, /,Cl;) are used to define the 
plasmon frequencies, %ujn- 

face, S, or bulk, B, determined by the rel- 
ative weight outside and inside the jellium 
perimeter, sketched in Fig. [3] (d). The po- 
tential and density of surface modes has a 
large contribution in the spill-out region out- 
side the jellium edge, and has lower energy 
than the bulk modes. The modes are labeled 
according to their angular quantum number 



m, accounting for the e^^^ angular symme- 
try, easily visualized for the surface modes 
as the number of standing waves along the 
wire perimeter. For the thinnest wire only 
modes up to m = 1 are sustained, whereas 
modes with larger m are found for increasing 
radii. This corresponds to a minimum plas- 
mon wavelength, Amin ^ 10 A that can exist 
for the structure, which is obtained by defin- 
ing an effective wave- vector for the modes, 
9eff = 9^ + Q'i- The parallel momentum trans- 
fer q\\ is fixed in the calculation, and q^ can 
be written: 



?± 



27r m 



(2) 



where A^ is the wavelength associated with 
the shape of modes perpendicular to the wire. 
When m is sufficiently large, or R small, the 
effective wave-vector becomes large enough 
for the mode to enter the regime of Landau 
damping (coupling to intra-band transitions 
with large q) , until the point where the mode 
can no longer be sustained. This is also re- 
flected in a mode strength, a, that decreases 
with m > 1 for the surface modes. 

For the bulk modes an additional quan- 
tum number, n, is required to account for 
radial dependence, where Bo^2 has an ex- 
tra node in the induced density compared to 
Bq^i. From a group-theoretical analysis of the 
possible symmetries within the point-group 
Z)oo/i5 the modes with m > are expected 
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(c) Radius 7.5 A 
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(d) Bulk and surface mode regions 



(e) Rotation symmetry and degeneracy 
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FIG. 3. Plasmon eigenmodes of Na jellium wires of radii: (a) 2 A, (b) 5 A, and (c) 7.5 A, 
corresponding to induced potential, (/>(r), and density, p(r). The strength, a, and energy, ?iuj^ is 
given for each mode. Modes with surface as well as bulk character are found, labelled S and B 
respectively, where the surface modes has a significant contribution in the spill-out region outside 
the jellium edge, see (d). The first subscript denotes the angular quantum number tti, describing 
the e^^^ rotational symmetry. The increase in m leads to a larger effective wave-vector for the 
plasmon, ^eff ? resulting in an increase in energy with m. Also the Landau damping increases with 
^eff , which is why the a larger number of modes can be sustained by wires with larger radius. As 
explained in (e), modes with tti > appear as degenerate pairs, rotated with T{/2m (plots labeled 
x2 when both are not shown). 



to appear in degenerate pairs. This is ex- identical modes with a relative rotation of 

actly found for the calculated modes, where 7r/2m, see Fig. [3] (e) 

pairwise identical eigenvalue- curves results in ^^ ^ave investigated the energy disper- 
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FIG. 4. Energy dispersion of plasmon eigenmodes for the quantum calculations (circles and 
squares) and the classical results for the surface modes (solid lines). The energy is plotted against: 
(a) Wire radius for fixed gy = O.IOA . (b) Wave- vector gy for fixed R — 5A. Circles and squares 
represent surface and bulk modes respectively, and the diameter/height of the points indicate the 
mode strength, a, obtained by fitting the eigenvalues to a single pole model, Eq. [T] The S^o mode, 
which has acoustic character, tends to zero energy for vanishing R or g, where there is a good 
agreement with the classical result. For large q the quantum model produces a positive dispersion 
for all modes, which is not captured by local models. For surface modes with tti > 0, there is an 
overall poor agreement due to effects from a larger g^gp and spill-out. 



sion of the modes with R and gy, shown 
in Fig. [4j where the result from our quan- 
tum calculation has been plotted together 
with the classical result obtained from solving 
Laplace's equation and using a local Drude 
model for the dielectric functioipS. The clas- 
sical model yields a dispersion relation for 
the surface modes that depends on a single 
dimensionless parameter, gi?, where q is the 
wave- vector parallel to the wire: 



^m = ^v ^^ ^L{Q^) Km{qR)- 



(3) 



Here /^ and K^ are the special Bessel's func- 
tions of order m, and uOp is the bulk plasmon 
frequency. The acoustic character of the S'o 
mode leads to an energy that goes to zero as 
a function of q\\R^ whereas the S'^>o modes 
tends towards the classical surface plasmon 
frequency ujpj\[2 for low q\\R. For large q\\R 
all modes tends towards ujpj^f^ for the clas- 
sical result. 

When comparing to our quantum results, 
a very good agreement is seen for the S^ 



mode at low q\\ , which breaks down for larger 
q\\. Since the Drude model assumes the elec- 
tronic response to be local in space, it cannot 
capture the dispersion with q. This explains 
the significant deviation for large gy, where 
the quantum result for the energy continue 
to increase. This g-dependence also explains 
why the energy of the surface modes increases 
with larger m (due to a larger Q'gff), opposed 
to the classical result where the energy does 
not exceed cjp/v^. However, Si and 82-, are 
seen to be red-shifted with respect to the clas- 
sical value in the limit Q' ^ 0, and for Si 
the redshift increases further with decreasing 
R. A redshift is expected for decreasing size 
of simple metal nanostructures due to elec- 
tron spill-out; as the surface to bulk ratio in- 
creases, the spill-out will lead to a decrease 
in the mean electron density, which results in 
a lower plasmon frequency. This redshift was 
also observed for the plasmon eigenmodes of 
thin slabs of sodiumP^, and for quantum cal- 
culation on small Na clusters^-*^ . This is op- 
posed to the behavior of silver nanostruc- 
tures, where coupling to single-particle tran- 
sitions originating from the d-states of silver 
leads to blueshift with decreasing siz d^^ ' ^^ l 
Thus, there are two competing factors affect- 
ing the energy of the S'^>o modes for decreas- 
ing R; the increase due to a larger perpen- 
dicular g, and the decrease due to spill out, 
which might explain the oscillatory behavior 




FIG. 5. The unit cell containing two jellium 
wires of distance d. We consider separations in 
the range to 20 A, and choose i? = 2.5 A for 
both wires, in order to study the case where only 
three modes exist for the individual wires. 

of the Si mode for small R. 

When considering the strength of the 
modes (represented by the diameter/height 
of the labels in Fig. lil), there is a clear ten- 
dency for weakening of modes for small R and 
large q\\. This is due to a smaller number of 
electronic bands crossing the Fermi level for 
the thin wires, which gives a smaller number 
of transitions contributing to the plasmon in 
this limit. The weakening with increasing q 
can be explained by increasing Landau damp- 
ing as discussed previously. 

After establishing an understanding of the 
quantum modes of single wires, we are ready 
to study the hybridization of modes between 
two parallel wires of small separation. The 
unit cell is sketched in Fig. [5| where the 
separation, rf, is varied from 20 A down to 
zero measured from the jellium edge, and the 



amount of vacuum separating the wires of In this regime of zero overlap between the 

neighboring super-cells are kept constant at electronic states, the modes can be classified 

20 A. The hybridization of plasmon modes according to conventional hybridization the- 

is not unlike the formation of molecular or- ory as even or odd combinations of the ini- 

bitals as bonding and antibonding combina- tial modes, labelled -|— I- and H — respectively, 

tions of atomic orbitala^DZZI^ where the new where the four S'l-types are further divided 

hybridized modes correspond to the even and into vertical, v, and horizontal, h, according 

odd combinations of isolated modes of the to the direction of polarization. As can be 

individual structures. At very small sepa- seen in Fig. [7| showing the induced potential 

rations however, this simplified description and density for 8 A, 2 A, and 1 A separa- 

can no longer hold since overlap of the elec- tions, this description is no longer valid for 

tronic states and eventual tunneling will alter small separations. 

the resonances. We have applied our quan- rp^^ ^^^^^^ dispersion of the modes as a 

tum approach to provide real-space visualiza- function of d is shown in Fig. gb), where 

tion of the hybridized modes in the transition ^^^ ^_ ^^^ ++ combinations split in energy 

from the non overlap to the overlap regime, for intermediate separations. A downshift is 

In order to hmit the number of hybridized found for modes with opposite charges facing 
modes we choose both wire radii equal to 2.5 each other, i.e. S^~ , 5J^++, and Sj+~ , which 
A, where only the Sq and the two degenerate can be explained by the electrostatic attrac- 
Si modes exists. Only modes of similar sym- tion between induced charges of opposite sign 
metry will couple to form hybridized modes; that leads to a lowering of the plasmon fre- 
for example, an 5*0 mode at one wire wiU be quency (and in reverse for the other modes), 
unable of inducing a S'l-like potential at a The coupling between the wires increases for 
neighboring wire, and no self-sustained oscil- decreasing separations which leads to a larger 
lation can occur. Therefore the modes only splitting, until the tunneling regime is en- 
couple to those of the same type {Si couple tered around 5 A. The shift is more distinct 
to the Si mode of same rotation), and six hy- for the S'o-types that also has the most long- 
bridized modes is predicted to exists for this ranged potential that explains the stronger 
structure. coupling and larger splitting of bonding and 

Six modes are found in the quantum calcu- antibonding modes. The potential from the 

lation as shown in Fig. [6|^a), with eigenmodes electric dipole of the 5'i modes decreases 

calculated for a 20 A separation of the wires, faster along y, in particular for the vertical 




FIG. 6. (a) Induced density of hybridized modes of two parallel wires with i? = 2.5 A at 20 A 
separation. The modes can be identified as hybridization between the So, Sf and Sf modes of 
the individual wires, where only modes of the same type couple due to symmetry restrictions, (b) 
The mode energies as a function of wire separation, d, with the mode strength indicated by the 
diameter of the points. The hybridization leads to a split in energy of the even, ++, and odd, -\ — 
mode-combinations, plotted together with the classical result (sohd curves) in (c) and (d). At 2 
A separation the Sq and S^ ^^ modes are quenched out due to tunneling and a charge-transfer 
mode Si^ appears. 



modes, which again shows a smaller splitting 
than the horizontal modes. 

In Fig. |6] (c) and (d) the quantum results 
have been plotted together with the results 
of a classical model. The later has been ob- 
tained by solving Laplace's equation for a sys- 
tem of two cylinders of equal radii R and di- 

electric constant e = 1 f . Following a 

general approach by A. Moradi^^, we have 
derived the expressions for hybridized Sq and 
S'l modes for this system, see Appendix. 

For separation > 15 A there is a good 
agreement between the classical and quan- 
tum results for the energy splitting. For de- 
creasing separations however, the quantum 



calculations produces a smaller splitting than 
the classical model, until all modes eventu- 
ally shows a positive dispersion with decreas- 
ing d. For the Sq modes the deviations from 
the classical model begins around d = 10 
A, whereas for the Si modes deviation is 
seen already at larger separation, approxi- 
mately at rf = 15 A. We explain the devi- 
ation by a beginning overlap of unoccupied 
electronic states of the two wires that are 
more long-ranged, and thus provides a vir- 
tual electronic tunneling. At 10 A separa- 
tion, bands from 1.7 eV above the Fermi level 
were found to be modified compared to 20 A 
separation, while noticeable changes in the 
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FIG. 7. The plasmon eigenmodes for two parallel nanowires with i? = 2.5 A for separations (a) 8 
A, (b) 2 A and (c) 1 A. Already at 8 A separation the potential, 0(r), and density, p(r), of the 
modes is slightly modified, the change being most dominant for the Si ^^ and Sq modes. At 2 A 
these mode are almost quenched out due to tunneling, and the induced density is shifted into the 
bulk. Below 2 A a new charge-transfer mode Sf^, corresponding to density oscillation between the 
wires, appears. For the other modes there is a continuous evolution of shape from the non-touching 
to the touching regime, as seen for the mode energies as well. 



occupied bands were only found below 4 A can be seen in Fig. [TJ already at 8 A sep- 
separation. aration some eigenmodes are modified com- 

pared to the large separation limit shown in 

Our approach for visualizing the eigen- r-. 

Fig. [6| (a). The modes are seen to be gradu- 

modes from a quantum calculation provides 

ally altered with decreasing separation, with 

a clear picture of how the plasmon modes , _ , ^v.^^ 

the Sq and Si modes being most seri- 

are altered due to these quantum eflFects. As 
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ously affected due to the facing of opposite shown how the hybridized plasmon modes of 
charges. The potential between the wires de- a two- wire system is modified due to quan- 
creases in amphtude going from 8 A to 2 A tum effects for decreasing separation. For 
separation, and the induced density is shifted separations above 15 A the modes were weU 
further into the bulk. A decrease in mode characterized according to conventional hy- 
strength (indicated by the marker diameter bridization theory, and a good correspon- 
in Fig. [6]) is also found for these modes upon dence was found with a classical model for 
approaching 2 A separation. At 1 A both the splitting of even and odd modes. This 
modes are quenched due to tunneling and are agreement was found to break down start- 
replaced by a dipolar charge-transfer mode, ing from 10 A to 15 A separation, where our 
here named /Sf^, and the induced density is quantum calculations produced a significant 
now located at the other side of the wires, smaller splitting. We attribute this to virtual 
away from the interface. That this mode is electronic tunneling due to overlap between 
different in nature than the Si ^^ and Sq~ unoccupied states, which in the range of 5- 
is also clear from Fig. [6](b), where it appears 15 A separation is strong enough to decrease 
as a separate line in the energy dispersion the induced density (and potential) in the in- 
just below 2 A, and increases in strength for terface, which results in a lower energy split- 
smaller separations. In this touching regime ting. Below 5 A the dispersion shows a qual- 
with large degree of tunneling, the structure itatively different behavior with an increase 
can be regarded as one (non-cylindrical) wire; in energy for all modes, indicating an larger 
however, for the other modes there is a con- degree of tunneling. However; only for sepa- 
tinuous transition from the non-touching to rations below 2 A did we find a true charge- 
the touching regime. The Sj^^ and Si ^~ transfer plasmon to appear. 



modes gradually build up induced density in 
the region between the wires, and Sj^~ ends 
up resembling a quadropolar (m=2) mode for 
the joint structure. 



Appendix: Classical hybridization 
model for a two-wire system 



We have derived the classical expression 
TV rONCTlJSTON ^^^ ^^^ energy of the hybridized modes of 

a system of two cylinders of equal radii i?. 

We have given a visual demonstration of where the general approach by A. MoradP^ 

the plasmon eigenmodes in nano-wires and simplifies to finding the zeros of the deter- 

12 



minant of the matrix: (M^ — /) , where the accounting for the sphtting of the So~ and 

elements of M are written: Sq~^ modes. Due to the smaU discrepancy 

u)'^ K (ad)! (qR) between the quantum results and the classi- 

cj2 — u)"^ Kn{qR) cal model for a single wire (due to dispersion 



where the m, n indices corresponds to angular '^i^^ Q ^^^ R), we have set ojq equal to the 
quantum numbers, and uj^ is the energy of result from the quantum calculation for Sq 
the Sm mode of a single wire given in Eq. |3| ^r a clearer comparison. For coupling of ^i 
For the couphng of two So modes both m,n= modes, m,n = ±1, and (M^ - /) is a 2 by 
0, and solving for omega yields two solutions: 2 matrix, which leads to a total of four solu- 
tions for uj: 



UJ. 






CUn 



l±Mqd) 



lojqR) 
KoiqR) 



(A.2) 



<^: 



{h,v},{±} 



= UJ. 



1{±, ^}^f^JK2{qdf + Koiqdf ± 2Ko{qd)K2{qd) 



(A.3) 



where uji has been set equal to the quantum 
result for Si. 
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